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Polarization  losses  within  electrodes  of  solid  oxide  fuel  cells  (SOFCs)  are  determined  both  by  material 
composition  and  microstructure.  Improvement  in  performance  can  be  supported  by  a  detailed  charac¬ 
terization  and  modeling  of  the  electrode  microstructure.  Focused  ion  beam  (FIB)  and  scanning  electron 
microscopy  (SEM)  combined  with  image  processing  have  already  proven  potential  for  the  reconstruction 
of  porous  electrodes.  In  this  contribution  the  serialized  reconstruction  procedure  of  a  high-performance, 
mixed  ionic-electronic  conducting  Lao.58Sr0.4Coo.2Feo.803_a  (LSCF)-cathode  will  be  illustrated  in  detail. 
Based  on  corrected  reconstruction  data  sets  and  by  the  evaluation  of  qualified  algorithms  discriminating 
between  porosity  and  electrode  material,  a  sensitivity  analysis  of  the  grayscale  threshold  value  on  the 
essential  parameters  (i)  surface  area,  (ii)  volume/porosity  fraction  and  (iii)  tortuosity  is  performed. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

To  accelerate  the  commercialization  of  solid  oxide  fuel  cells 
(SOFCs),  it  is  essential  to  improve  the  electrode  performance  over 
the  entire  temperature  range  between  600  and  900  °C.  Mixed  ionic- 
electronic  conductors  (MIECs)  such  as  La0.58Sro.4Coo.2Feo.803_5 
(LSCF)  are  qualified  as  high-performance  cathodes  at  interme¬ 
diate  and  low  operating  temperatures.  Polarization  losses  are 
determined  both  by  material  composition  and  microstructure. 
Improvement  in  performance  can  be  supported  by  a  detailed 
characterization  and  modeling  of  the  electrode  microstructure. 
High-quality,  detailed  microstructure-data  are  a  prerequisite, 
which  can  be  obtained  by  focused  ion  beam  (FIB)  and  scanning 
electron  microscopy  (SEM)  combined  with  image  processing  [1-6]. 
The  FIB  procedure  applies  Ga  ions  to  mill  thin  sections  from  the 
exposed  electrode  surface,  which  is  imaged  at  the  same  time  by 
SEM.  Sequential  milling  and  imaging  yield  a  serial  set  of  consec¬ 
utive  images  of  the  electrode.  By  stacking  the  2D  images  in  a  3D 
space  (cf.  Fig.  1 )  and  aligning  the  pixels  of  the  images  in  the  slicing 
direction  (space  between  two  images),  a  3D  reconstruction  of  the 
microstructure,  consisting  of  voxels  (volumetric  pixels),  is  derived 
[6]. 
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This  3D  reconstruction  data  set  gives  a  valuable  basis  for  the 
determination  of  essential  microstructure  parameters:  (i)  surface 
area,  (ii)  volume/porosity  fraction  and  (iii)  tortuosity.  This  com¬ 
bined  FIB/SEM  and  reconstruction  approach  was  reported  recently 
for  both  SOFC  anodes  and  cathodes  by  several  groups  [1-5]. 
However,  the  reconstructed  volume  elements  were  kept  rather 
small  and  hardly  representative  of  the  entire  electrode  structures. 
Furthermore,  aspects  concerning  the  accuracy  of  the  alignment 
procedure  were  not  given  in  detail,  but  are  essential  for  the  quality 
of  the  3D  reconstruction  data  sets,  as  discussed  in  [7].  The  unam¬ 
biguous  segmentation  of  the  grayscale  frequency  distribution  into 
pores  and  electrode  material,  to  our  understanding,  is  the  most 
important  step  in  image  processing.  Therefore,  we  will  present  (1 ) 
an  analysis  of  qualified  algorithms  able  to  automatically  identify  a 
“true”  threshold  value  and  (2)  discuss  the  variance  of  the  resulting 
microstructure  parameters. 

2.  Materials  and  methods 

2.1.  Sample  preparation 

The  FIB/SEM  procedure  was  applied  to  a  high-performance 
LSCF  cathode,  which  was  screen-printed  and  sintered  onto  a 
Ce0.8Gdo.202_5  (CGO)  interlayer  and  a  thin  film  Y0.i6Zr0.84O2_5  elec¬ 
trolyte  (8YSZ).  The  entire  anode  supported  cell  (ASC)  was  developed 
at  the  Research  Center  Jiilich,  details  of  the  manufacturing  proce¬ 
dure  are  given  in  [8].  Moreover,  the  electrical  performance  of  the 


0378-7753 /$  -  see  front  matter  ©  2010  Elsevier  B.V.  All  rights  reserved, 
doi:  1 0.1 01 6/j.jpowsour.201 0.1 0.006 


J.Joos  et  al.  /  Journal  of  Power  Sources  196(2011)  7302-7307 


7303 


al 


Fig.  1.  SEM-images  of  a  porous  Lao.58Sr0.4Coo.2Feo.803_5  (LSCF)  cathode  (al  and  a2):  the  usable  sector  (bl  and  b2)  of  the  entire  image  shifts  upwards  with  increasing  penetration 
depth  in  slice  direction  (z-direction).  By  stacking  and  aligning  the  2D  images  in  3D  space  a  detailed  3D  reconstruction  of  the  microstructure  is  derived  (d)  [6]. 


ASC  was  evaluated  in  the  past  in  great  detail  via  electrochemical 
impedance  spectroscopy  [9,10]. 

The  sample  preparation  was  presented  in  detail  in  a  previous 
paper  [6],  nevertheless  the  most  important  issues  will  be  repeated 
here.  The  ASC  was  first  cut  into  smaller  pieces  of  approximately 
5  mm  x  2  mm.  Then  the  pores  were  vacuum  infiltrated  with  a  two 
component  resin,  which  strengthens  the  porous  microstructure 
during  FIB/SEM  treatment  and  allows  for  a  better  planar  sectioning. 
But  most  important  it  supports  the  identification  of  the  pore  phase 
during  image  processing  (as  the  resin  within  the  pores  appears  as 
dark  black,  which  improves  the  contrast  with  the  light  gray  elec¬ 
trode  material,  see  SEM  images  in  Fig.  1 ).  Hence  infiltration  is  an 
essential  precondition  to  allow  an  automated  segmentation  on  the 
consecutive  images,  which  will  be  discussed  in  detail  later. 

The  edges  of  the  sample  were  polished  to  obtain  planar  surfaces 
and  sharp  edges.  This  allows  starting  the  process  of  FIB  section¬ 
ing  right  from  the  edge  of  the  sample.  A  gold  coating  (~2.5  p,m 
thickness)  was  sputtered  on  top  of  the  sample  to  protect  the 
microstructure  and  help  to  reduce  streaking  [2]. 

2.2.  FIB/SEM  processing 

For  the  sequential  process  of  FIB  milling  and  SEM  imaging  a 
ZEISS  1540XB  CrossBeam®  (Carl  Zeiss  AG,  Oberkochen,  Germany) 
was  used.  The  two  beams  have  a  coincident  angle  of  54°,  as  shown 
schematically  in  Fig.  2.  It  is  possible  to  run  the  sequential  process 
automated,  but  before  starting  the  process,  it  is  important  to  find 
the  right  settings.  Therefore  we  made  some  preliminary  calcula¬ 
tions  concerning  the  accuracy  of  the  data  obtained  with  FIB/SEM 


techniques  [6].  An  important  factor  is  the  resolution  or  rather  the 
number  of  pixels  (respectively  voxels)  per  particle-diameter  that  is 
necessary  in  the  SEM  images.  The  results  [6]  have  shown  that  a  res¬ 
olution  of  at  least  20-30  voxels  per  particle-diameter  is  required 
for  the  convergence  of  parameters  like  surface  area  or  porosity. 
To  ensure  a  sufficient  resolution  in  all  three  directions,  this  holds 
not  only  for  the  resolution  of  the  SEM  images  (resolution  in  x-  and 
y-direction)  but  also  for  the  slicing  resolution  (distance  between 
two  images,  i.e.,  resolution  in  the  z-direction).  However,  in  most 
cases  there  is  a  difference  between  the  resolution  of  the  SEM  images 
and  the  slicing  distance  [1-4].  This  can  be  readjusted  by  applying  a 
re-sample  step  on  the  FIB/SEM  data. 

2.3.  3D  reconstruction 

An  important  aspect  is  the  estimation  of  the  minimal  volume 
size  that  has  to  be  reconstructed  for  reliable  and  representative 
results.  This  volume  size  is  generally  named  representative  volume 
element  (RVE)  or  stochastic  equivalent  representative  volume  ele¬ 
ment  (SERVE),  and  this  aspect  was  discussed  in  detail  in  [6,11,12]. 
Here  a  SERVE  is  employed  but  we  use  the  term  RVE  in  accordance 
with  other  publications  [7,11].  Based  on  a  coarse  3D  FEM-model, 
described  in  [6,13,14],  we  calculated  the  area  specific  resistance  of 
the  cathode  ( ASRcat )  as  a  figure  of  merit  for  cathode  performance. 
As  a  result  the  ASRcat  drops  with  increasing  number  of  electrode 
particles  N  (hence  with  increasing  contact  area  between  cathode 
and  electrolyte)  and  reaches  an  asymptotic  limit.  Logically,  this 
limit  depends  on  the  porosity  of  the  electrode  structure  consid¬ 
ered.  In  the  case  presented,  as  the  porosity  is  in  between  30  and 


Fig.  2.  (a)  Schematic  design  of  a  FIB/SEM  setup  from  the  side  view  [3].  (b)  Prepared  sample  and  (c)  the  sample  with  the  observed  region  of  interest  after  FIB/SEM  treatment. 
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Fig.  3.  Histogram  of  the  grayscale  frequency  distribution  of  the  reconstruction  data: 
thin  blue  lines  show  the  individual  distributions  of  grayscales  for  each  of  the  213 
images  used  for  the  reconstruction,  while  lines  in  red,  green  and  ochre  indicate  the 
average  values  of  3  RVEs.  The  appropriate  threshold  value  for  the  segmentation 
into  pores  and  electrode  material  (LSCF)  is  located  in  the  “valley”  between  the  two 
histogram  peaks. 

50%,  N  should  be  >6-7  in  one  dimension  [6].  We  conclude  that 
the  volume  element  must  cover  a  base  area  of  at  least  7x7  parti¬ 
cles  [6,13].  The  reconstructed  cathode  structure  shown  in  Fig.  Id  is 
30.7  pan  (height)  x  5.25  pan  x  15.75  pan,  and  was  assembled  from 
almost  seven  hundred  consecutive  images. 

The  SEM  images  taken  during  FIB  milling  are  usually  comprised 
of  a  larger  area  than  required  for  reconstruction.  During  FIB  milling, 
the  entire  image  shifts  upwards  with  increasing  penetration  depth 
in  the  slicing  direction  (z-direction).  Therefore,  the  usable  sector  for 
reconstruction  has  to  be  cut  out  of  the  entire  image.  Fig.  1  shows  the 
SEM-images  of  the  investigated  porous  LSCF-cathode  (Fig.  1,  Parts 
al  and  a2)  as  well  as  the  sector  used  for  reconstruction  (Fig.  1 ,  Parts 
bl  and  b2).  The  cut-out  procedure  was  supported  by  the  software 
MATLAB  (The  MathWorks,  Natik,  MA,  USA).  By  stacking  all  seven 
hundred  2D  images  in  a  3D  space  and  expanding  the  pixels  in  the 
third  dimension  (space  between  two  images),  the  pixels  become 
voxels  and  a  3D  reconstruction  is  obtained.  Thereupon  the  precise 
alignment  of  consecutive  images  is  of  outstanding  importance  for 
the  resulting  3D  reconstruction  data  set.  So  far,  this  issue  was  not 
extensively  discussed  in  the  literature.  In  our  case,  even  small  mis¬ 
alignment  leads  to  defective  information  due  to  the  high  resolution 
of  the  images  (pixel  size  of  ~35  nm,  image  distance  of  ~25  nm),  for 
more  details  we  refer  to  [7]. 

3.  Results  and  discussion 

3.1.  Threshold  grayscale  value 

One  single  resulting  3D  reconstruction  data  set,  denominated 
as  a  RVE,  contains  213  individual  images.  Consequently  we  could 
generate  three  non-overlapping  RVEs  out  of  seven  hundred  images. 
Every  single  image  consists  of  different  grayscale  values,  discrim¬ 
inated  from  a  value  of  0  (black)  to  a  value  of  255  (white).  The 
resulting  grayscale  value  frequency  distribution  of  all  213  images 
as  well  as  the  average  values  of  the  3  RVEs  are  shown  in  the  his¬ 
togram  (Fig.  3).  At  a  glance,  the  most  significant  information  is 
extracted  from  the  two  maxima,  clearly  assigned  to  (i)  pores  at  a 
mean  grayscale  value  of  ~40  and  (ii)  electrode  material  at  ~144. 

More  precisely,  the  entire  grayscale  image,  voxel  by  voxel,  has 
to  be  unambiguously  assigned  either  to  a  pore  or  to  electrode  mate¬ 
rial.  This  so-called  segmentation,  to  our  understanding,  is  the  most 
important  step  in  image  processing  and  potentially  a  large  source  of 


error.  For  segmentation,  the  literature  proposes  different  methods 
and  it  is  even  today  a  field  of  ongoing  research.  For  the  segmenta¬ 
tion  of  pixel-  and  voxel-based  data,  the  most  common  and  capable 
method  of  image  segmentation  is  thresholding  [15].  The  appropri¬ 
ate  threshold  value  for  the  reconstruction  is  located  somewhere  in 
the  “valley”  between  the  two  histogram  peaks  (cf.  Fig.  3). 

But  despite  its  apparent  simplicity,  thresholding  is  a  very  power¬ 
ful  method  but  also  challenging  [15].  First  and  foremost  it  has  to  be 
clarified  that  segmentation  with  a  threshold  value  is  a  constructive 
approach,  followed  by  the  allocation  of  the  “true”  threshold  value. 
To  check  if  an  automated  segmentation  with  a  threshold  value  is 
expedient,  a  shape  analysis  of  the  histogram  (Fig.  3)  can  be  helpful. 

A  strong  indication  that  segmentation  via  thresholding  with 
one  threshold  value  is  admissible  is  a  bimodal  distribution  of  the 
grayscale  values,  which  is  obviously  fulfilled  here.  Poor  image  con¬ 
trast  leads  to  a  histogram  where  peaks  tend  to  overlap,  which 
makes  it  difficult  to  separate  pore  from  electrode  material.  Even 
more  serious  are  image  shadowing  artifacts  [15],  which  lead  to 
images  that  appear  lighter  at  one  side  and  darker  at  the  other.  Fur¬ 
thermore,  because  of  the  need  for  a  reasonable  number  of  images 
(here:  one  RVE  consists  of  213  images),  all  of  them  should  have 
the  same  contrast.  Additionally  the  absence  of  a  brightness  gradi¬ 
ent  between  the  consecutive  images  is  a  precondition.  This  can  be 
controlled  by  comparing  the  individual  histograms  of  the  subse¬ 
quent  images.  If  the  histogram  peaks  of  the  different  images  are 
not  shifted  against  each  other,  it  can  be  concluded  that  there  is  no 
brightness  gradient  and  no  image  shadowing. 

It  can  be  recognized  from  the  grayscale  frequency  distribution 
of  213  individual  images  (cf.  Fig.  3)  that  the  distribution  of  the 
two  peaks  is  well  aligned.  As  a  consequence,  brightness  gradient 
or  shadowing  is  of  no  concern  and  segmentation  with  a  uniform 
threshold  value  is  permitted.  (In  our  case,  image  shadowing  was 
avoided  by  an  advantageous  sample  preparation  [6],  alternatively 
the  lift-out  technique  can  be  applied  [2].) 

As  the  preconditions  for  segmentation  by  thresholding  are  now 
fulfilled,  we  can  concentrate  on  the  question:  which  grayscale  value 
in  the  “valley”  represents  the  appropriate  threshold?  Locating  the 
“true”  threshold  value  is  the  most  important  task  in  image  pro¬ 
cessing.  Fig.  4  shows  two  examples  for  a  poorly  chosen  threshold: 
if  the  corresponding  grayscale  value  is  much  too  low  (with  a  value 
60  in  Fig.  4a),  a  significant  part  of  the  pore  volume  is  attributed 
to  electrode  material,  and  if  the  corresponding  grayscale  value  is 
chosen  too  high  (with  a  value  of  120  in  Fig.  4b),  electrode  material 
may  be  attributed  to  pore  volume.  It  is  obvious  that  the  resulting 
porosity /volume  fraction  would  be  inaccurate  for  both  cases. 

The  “true”  threshold  can  be  roughly  allocated  by  interpreting 
the  histogram  shape.  Given  a  bimodal  distribution,  the  minimum 
or  “valley-floor”  between  the  peaks  is  the  obvious  choice  for  the 
threshold.  In  our  case,  where  the  two  peaks  are  not  clearly  resolv¬ 
able,  several  algorithms  were  developed  during  the  last  years.  Most 
of  them  are  based  on  statistical  interpretation  techniques  and  auto¬ 
matically  locate  the  “true”  threshold. 

3.2.  Algorithm  for  threshold  determination 

As  in  our  histogram  the  peaks  are  not  clearly  resolved,  respec¬ 
tively  no  clear  minimum  is  observed  (cf.  Fig.  3),  we  evaluated 
several  methods  which  automatically  identify  the  “true”  thresh¬ 
old  value.  A  widely  used  method  is  the  Ptile-  or  Median-algorithm 
[16].  This  algorithm  works  well  if  an  a-priori  knowledge  of  the  vol¬ 
ume  fraction  is  available.  This  method  is  obviously  inadequate,  as 
the  determination  of  the  volume  fractions  is  one  of  the  goals  of  this 
work.  A  very  intuitive  method  for  finding  the  appropriate  thresh¬ 
old  is  the  Mean-algorithm  [17],  which  calculates  the  mean  value  of 
the  grayscale  values.  This  algorithm,  when  applied  on  our  data  set, 
determines  a  threshold  value  of  100. 
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Fig.  4.  Examples  for  a  poorly  chosen  threshold:  (a)  if  the  threshold  is  chosen  too  low,  parts  of  pores  are  attributed  wrongly  to  electrode  material  (dark  green  regions),  while 
(b)  if  the  threshold  is  chosen  too  high,  parts  of  electrode  material  are  attributed  wrongly  to  pores  (light  gray  regions).  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


A  more  advanced  and  one  of  the  most  qualified  methods  for 
image  data  consisting  of  two  phases  (pore  and  electrode  material) 
is  Otsu’s  method  [  1 8  ].  The  two  phases  are  considered  as  two  clusters 
Ci  and  C2  of  grayscale  values,  which  are  defined  as 


255 


Ci  =  y^p(i)andC2  =  Y  p(i). 


i=0 


i=T+ 1 


(1) 


The  algorithm  sets  both  clusters  as  tightly  as  possible.  While 
adjusting  the  threshold  value  the  spread  of  one  cluster  increases, 
whereas  the  other  decreases.  The  target  is  then  to  minimize  the 
combined  spreads  or  equivalently  maximize  the  variance  between 
the  clusters.  First,  the  variances  of  the  two  clusters  of  intensity  val¬ 
ues,  separated  by  an  intensity  threshold  T,  are  calculated  by  using 
the  mean  values  ii\  and  fi2  of  their  grayscales: 


T  255 

a\  =  YU  -  fii)2p(i)anda%  =  y~)(i-  /x2)2p(i).  (2) 

i=o  i=r+i 


Second,  the  variance  between  the  clusters  is  calculated  by  using  the 
mean  value  of  all  grayscale  values  /x: 

&Bj  =  Cl (/Xi  -  p.)2  +  C2(p.2  —  AO2-  (3) 

Now,  the  function 

-*  max  (4) 

ffT 

is  calculated  for  every  intensity  i  =  0,. . .,/  -1,  where  <7j  =  C\o\  + 
C2crj  is  the  total  variance.  The  intensity  which  maximizes  this  func¬ 
tion  (4)  is  said  to  be  the  optimal  threshold. 

We  applied  the  Otsu-algorithm  on  our  image  data  for  all  three 
directions  (x-,  y-  and  z-plane),  and  came  up  with  the  “true” 
threshold  value  of  97  for  all  three  directions.  The  threshold  val¬ 
ues  resulting  from  the  Otsu-algorithm  (mean  threshold  value  97) 
and  the  Mean-algorithm  (mean  threshold  value  100)  are  rather 
close  to  each  other,  thus  the  calculated  porosity  differs  only  by 
1%  (absolute).  Nevertheless,  we  favor  the  Otsu-algorithm  because 
this  algorithm  makes  use  of  more  information  by  calculating  the 
variance  of  the  two  phases.  Interestingly  enough,  we  extracted  the 
same  threshold  value  of  97  also  by  “intuitive”  inspection  of  the 
histogram.  In  the  end,  we  performed  the  segmentation  procedure 
automatically  with  a  threshold  value  of  97.  Thereby  we  fixed  a 
porosity  fraction  of  ~47%  (47.27,  see  Table  1)  for  the  porous  LSCF 
cathode,  which  is  on  the  upper  boundary  of  what  we  subjectively 
expected  beforehand.  A  detailed  description  of  the  quantification 
of  microstructure  parameters,  as  well  as  a  sensitivity  analysis  on 
the  threshold  value  will  be  given  in  the  next  section. 


3.3.  Quantification  of  microstructure  parameters 

From  the  segmented  three-dimensional  data,  the  quantifica¬ 
tion  of  the  essential  parameters  (i)  volume/porosity  fraction,  (ii) 
electrode  surface  area  and  (iii)  tortuosity  of  pores  and  electrode 
material  was  realized.  However,  it  seems  that  there  is  no  estab¬ 
lished  method  to  evaluate  these  parameters  from  FIB/SEM  data  [3]. 
But  irrespective  of  the  method  by  which  the  parameters  were  cal¬ 
culated,  it  turned  out  that  the  results  depend  predominantly  on  the 
chosen  threshold.  Fig.  5  shows  the  course  of  electrode  surface  area 
a  and  porosity  fraction  e  in  dependence  of  threshold  values  varied 
from  60  to  100. 

As  expected,  porosity  is  highly  sensitive  to  the  chosen  thresh¬ 
old  value.  As  discussed  in  the  previous  Section,  with  increasing 
threshold  value,  more  and  more  voxels  are  attributed  to  pores,  thus 
leading  to  an  increase  of  the  porosity  fraction.  With  a  threshold  of 
60  the  porosity  is  clearly  below  35%  for  all  three  RVEs,  whereas  with 
a  threshold  of  100  the  porosity  increased  to  50%.  We  conclude  that 
it  is  essential  to  establish  the  “true”  threshold  value  for  the  calcula¬ 
tion  of  the  porosity /volume  fraction.  In  contrast,  the  surface  area  is 
nearly  insensitive  to  a  variation  in  threshold  value,  as  it  varies  only 


Table  1 

Grayscale  threshold  values  calculated  from  image  data  by  (a)  the  Otsu-algorithm 
and  (b)  the  Mean-algorithm  and  the  corresponding  figures  for  porosity  fraction  of 
one  RVE. 


Algorithm 

Threshold  value 

Porosity  fraction,  % 

Mean-algorithm 

100 

48.50 

Otsu-algorithm 

97 

47.27 

threshold  value 

Fig.  5.  Sensitivity  analysis  of  the  variance  in  threshold  values  on  the  microstructure 
parameters:  (i)  porosity  fraction  and  (ii)  electrode  surface  area,  calculated  with  3 
RVEs. 
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Fig.  6.  Left:  schematic  design  of  the  finite  element  method  model  which  was  applied  for  the  calculation  of  the  tortuosity.  Right:  example  of  the  potential  distribution  0  inside 
the  reconstructed  microstructure. 


in  between  6.2  (|xm_1 )  and  6.5  (|xm_1 )  which  is  not  surprising  for 
two-phase  materials. 

Before  calculating  the  tortuosity  of  the  electrode  structure, 
the  3D  reconstruction  data  were  re-sampled  to  obtain  cubic 
voxels  (35  nm  x  35  nm  x  35  nm).  Furthermore,  a  finite  element 
method  (FEM)  model  was  implemented,  which  has  also  been 
described  in  [7].  A  schematic  design  of  the  model  is  given  in 
Fig.  6.  The  tortuosity  is  determined  in  y-direction,  from  the 
electrolyte/cathode-interface  to  the  current  collector.  On  top  and 
underneath  the  reconstructed  volume,  two  thin  layers  with  a  high 
electrical  conductivity  crdist  are  applied  to  distribute  and  collect  the 
current  to  the  porous  electrode  structure  (cf.  Fig.  6).  A  potential 
difference  r\modei  between  these  distributing  layers  is  predefined  as 
1 .0  V  by  applying  a  potential  0  =  1 .0  V  on  the  upper  layer  and  0  =  0  V 
on  the  bottom  layer.  Inside  the  solid  phase,  the  transport  equation: 

V-(-^u/fcV0)  =  O  (5) 

is  solved  to  determine  the  current  /  that  circulates  through  the 
structure.  Thereby  the  intrinsic  conductivity  of  the  material  crbuik 
is  predefined.  At  the  interface  between  the  two  phases,  material 
and  pore,  a  no-flux  boundary  condition  was  applied.  Afterwards 
the  resistance  R  of  the  structure  is  calculated  by  using  Ohm’s  law. 
Hence  the  effective  conductivity  is  given  as 

J  Ny  .  . 

where  Nx,  Ny  and  Nz  are  the  number  of  voxels  in  x-,  y-  and  z- 
direction  and  lv  is  the  edge  length  of  the  cubic  voxels.  Finally  the 
tortuosity  of  the  material  (rmiec)  can  be  calculated  by  the  formula: 

W  =  (l-e)^,  (7) 

aeff 

where  s  is  the  previously  determined  porosity  fraction. 

To  determine  the  tortuosity  of  the  pores  (r poreX  an  artificial 
intrinsic  conductivity  obvdk  is  assigned  to  the  pores.  Therewith 
equation  (5)  is  solved  in  the  domain  of  the  pores  and  analogously 


to  the  previously  described  method,  the  tortuosity  of  the  pores  is 
given  by 


Tpore 


G bulk 
°eff  ’ 


(8) 


Based  on  this  procedure,  the  tortuosity  of  the  electrode  mate¬ 
rial  and  pores  was  calculated  for  the  threshold  values  of  70,  85 
and  97  and  then  averaged  from  3  RVEs.  Given  a  base  area  of 
5.25  [xm  x  5.25  |xm,  which  is  equivalent  to  152  x  152  voxels,  each 
of  the  three  RVEs  contains  ~1 8.5  x  1 06  voxels  for  a  height  of  28  |xm 
(which  is  the  layer  thickness  of  the  analyzed  LSCF  cathode).  The  fig¬ 
ure  for  tortuosity  had  been  solved  by  parallel  calculation  on  a  high 
performance  computer,  and  is  summarized  in  Table  2. 

As  can  be  seen,  the  tortuosity  changes  from  1.66  for  a  threshold 
of  70-2.06  for  a  threshold  of  97.  One  should  keep  in  mind,  that 
the  threshold  of  97  (tortuosity  2.06)  was  identified  by  the  Otsu- 
algorithm.  Clearly,  the  tortuosity  is  almost  as  sensitive  to  a  variation 
in  threshold  value  as  the  porosity.  In  other  words,  with  increasing 
pore  or  electrode  material  fraction,  the  tortuosity  of  the  contrary 
phase  (electrode  material  or  pore)  rises. 

Finally  we  calculated  the  microstructure  parameters  poros¬ 
ity,  surface  area  and  tortuosity,  for  a  (constant)  base  area  of 
5.25  [xm  x  5.25  |xm,  but  for  stepwise  increasing  height.  In  other 
words:  the  reconstructed  volume  increases  in  size,  starting  at 
1.75  [xm  and  ending  at  28  |xm,  which  is  almost  the  entire  layer 
thickness  of  the  analyzed  LSCF  cathode.  This  allows  us  to  study 
a  possible  influence  of  the  actual  reconstruction  volume  on 
the  microstructure  parameters.  The  results  are  summarized  in 
Table  3. 

Consequently,  for  the  porous  LSCF  structure  investigated  in 
this  paper,  the  volume  considered  should  take  at  least  half  of  the 
layer  thickness  into  account.  At  a  thickness  or  height  less  than 
h  <  14  |xm  we  do  not  obtain  a  representative  volume  element,  as  the 
microstructure  parameters  quite  obviously  then  depend  on  the  ele¬ 
ment  height  h.  This  dependency  disappears  for  larger  volumes,  as 
expected  for  homogeneous  microstructures.  In  any  case,  the  calcu¬ 


lable  2 

Average  tortuosity  of  the  electrode  material  and  pores  calculated  for  different  threshold  values  (average  of  3  RVE,  each  RVE  has  a  reconstructed  volume  size  of 
5.25  |jim  x  5.25  |jim  x  28  pirn,  or  ~18.5  x  106  voxels). 


Threshold  value  (TH)  (grayscale  frequency  distribution  spans  from  0  to  255) 

TH  70 

TH  85 

TH  97  (Otsu) 

Tortuosity  of  electrode  material  (LSCF)  Tmiec 

1.66 

1.89 

2.05 

Electrode  material  (LSCF)  fraction  (1  -e),% 

62.26 

55.91 

52.40 

Tortuosity  of  pore  r pore 

2.55 

2.11 

1.94 

Porosity  fraction  e,  % 

37.74 

44.09 

47.60 

J.Joos  et  al.  /  Journal  of  Power  Sources  196(2011)  7302-7307 


7307 


Table  3 

Microstructure  parameters  porosity,  surface  area  and  tortuosity,  calculated  for  a  (constant)  base  area  of  5.25  |xm  x  5.25  p,m  but  for  increasing  height  (reconstructed  volume 
has  increasing  size,  as  28  p,m  is  almost  the  entire  layer  thickness  of  the  analyzed  LSCF  cathode)  (threshold  value  of  97). 


Height  h,  (Jim 

1.75 

3.5 

7.0 

10.5 

14.0 

24.5 

28.0 

Porosity  fraction  e,  % 

48.91 

47.87 

47.82 

47.65 

47.81 

47.34 

47.27 

Surface  area  a,  |xm_1 

6.597 

6.674 

6.656 

6.502 

6.388 

6.276 

6.308 

Tortuosity  of  MIEC  xmiec 

1.954 

1.939 

1.950 

1.969 

1.989 

1.989 

1.985 

Tortuosity  of  pore  zpore 

2.039 

1.983 

1.975 

1.979 

1.959 

1.964 

1.967 

lations  must  be  solved  by  parallel  calculation  on  a  high  performance 
computer. 

4.  Conclusions 

The  microstructure  of  porous  electrodes  can  be  assessed  only 
if  accurate  data  on  grain  size,  porosity,  tortuosity,  etc.  are  avail¬ 
able.  For  this  purpose,  a  serialized  reconstruction  procedure  for 
FIB/SEM  derived  3D  microstructure  data  of  a  high-performance, 
mixed  ionic-electronic  conducting  Lao.58Sr0.4Coo.2Feo.803_«$  (LSCF)- 
cathode,  was  developed. 

During  the  entire  reconstruction  process  numerous  sources  of 
error  can  occur  ( see  [  6,7  ] )  but  in  this  contribution  we  focused  on  the 
section  of  image  processing.  Flerein,  the  unambiguous  segmenta¬ 
tion  of  the  grayscale  frequency  distribution  into  pores  and  electrode 
material,  to  our  understanding,  is  the  most  important  step.  There¬ 
fore,  an  analysis  of  two  qualified  algorithms  (Otsu,  Mean)  capable 
of  automatically  identifying  a  “true”  grayscale  threshold  value 
was  performed.  In  case  of  a  bimodal  (here:  pores  and  electrode 
material)  frequency  distribution,  Otsu’s  method  was  marginally 
preferred.  Furthermore  the  variance  of  the  parameters:  (i)  surface 
area,  (ii)  volume/porosity  fraction  and  (iii)  tortuosity  depending 
on  the  threshold  value  were  calculated.  Whereas  porosity  as  well 
as  tortuosity  varies  strongly  with  the  threshold  value,  the  sur¬ 
face  area  is  almost  constant  between  porosities  of  35  and  50%. 
For  the  given  LSCF  electrode  structure,  a  “true”  threshold  of  97 
(among  a  grayscale  variation  from  0  (black:  pore)  to  255  (white: 
electrode  material)  was  identified.  The  microstructure  parameters 
were  identified  as  (i)  surface  area  6.3  (pan-1 ),  (ii)  porosity  fraction 
47.3%  and  (iii)  tortuosity  of  pore  1.97. 

Part  of  the  results  could  only  be  extracted  from  a  newly  devel¬ 
oped  3D  FEM-model,  which  can  cope  with  a  large  number  of  voxels 
(107  to  108).  It  contains  a  solver  based  on  high  performance  com¬ 
puting  techniques,  and  a  3D  stochastic  geometry  generator.  Further 
details  on  this  3D  numerical  model  will  be  the  subject  of  a  forth¬ 
coming  article. 
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